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Abstract 



The growth by condensation of small water droplets in a three-dimensional ho- 
mogeneous isotropic turbulent flow is considered. Within a simple model of ad- 
vection and condensation, the dynamics and growth of millions of droplets are 
integrated. A droplet- size spectra broadening is obtained and it is shown to in- 
crease with the Reynolds number of turbulence, by means of two series of direct 
numerical simulations at increasing resolution. This is a key point towards a proper 
evaluation of the effects of turbulence for condensation in warm clouds, where the 
Reynolds numbers typically achieve huge values. The obtained droplet-size spec- 
tra broadening as a function of the Reynolds number is shown to be consistent with 
dimensional arguments. A generalization of this expectation to Reynolds numbers 
not accessible by DNS is proposed, yielding upper and lower bounds to the actual 
size-spectra broadening. A further DNS matching the large scales of the system 
suggests consistency of the picture drawn, while additional effort is needed to eval- 
uate the impact of this effect for condensation in more realistic cloud conditions. 
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1. Introduction 



The growth of microdrop lets by condensation is a long-stand i ng p roblem of cloud physics 
( Pruppacher and Klettl d igQVh). meteo rology (see e.g. llPCC-WGll jjoOlb ). medicine dMartonen 
(12000)) and engineering (iZhao et al.l (|l999h ). A fundamental understanding of key issues such 
as the turbulent mixing inside clouds, or the interaction of turbulence with microphysics is im- 
portant for a variety of applications (the parameterization of small scales in large scales models, 
the analysis of radiative transfer through clouds, the accurate prediction of the initiation of 
precipitation). The peculiar features of turbulence, influencing the motion on a wide range of 
space-time scales, can alter significantly the condensation process, providing a strongly fluctuat- 
ing and intermittent moist environment. This is a well-known concept in engineering, where tur- 
bulence is considered as a key i ngredient for the o ptimization of fuel-air mixing and of the rate 
of fuel vaporization in engines (IZhao et al.l (|l999|) ). Similar ideas have a long history in clouc 



physic s : turbulence plays a role both i n collision processe s (fSaffman and Turneil d 19561): Shaw 



(2003); S undaram and Collins (1 1997h : iGhosh et al.l (120051) : iFalkovich et all (l2002l) : IWang et al. 
(12006) : iWilkinson et al.l (|2006|) ) and in condensation. In this work we will focus particularly on 
the latter process. 

Condensation is a fundamental process for the early stages of cloud evolution. It is the 
only mechanism providing a growth of cloud droplets immediately after their formation, when 
they are few micrometers in size. The initiation of collisions and coalescence occurs when a few 
droplets become large enough to fall - a radius of 20 jim is commonly considered as a threshold. 
For these droplets, collisions are much more efficient than condensation as a growth mechanism. 
Indeed, after nucleation, droplet radius grows roughly one order of magnitude by condensation. 
In a timescale comparable to that of condensation, collisions produce raindrops, bridging a gap 
of about two orders of magnitude in size. However, the high efficiency of the latter is strongly 
influenced by the general features of the former. In particular, gravitational collisions are highly 
effective when the previous condensation process produces a population of droplets spanning 
a large variety of sizes. Unifo r m condensational growth l eads to narrow droplet-size spectra 
instead (see e.g. ISquiresI (ll952|) : lLee and Pruppacher! (| 19771) ). Indeed, provided that all droplets 
grow in similar ambient conditions, small droplets grow faster than large ones and thus all 
droplets finally tend to converge to the same size. This yields the long-standing problem of the 
bottleneck between condensation and collision-coalescence. 

The presence of broad droplet-size distribu tions has been detected via in situ measurements 
in warm clouds under very different conditions, Wameil (ii969). Still, this observation eludes full 
theoretical understanding despite the number of different approaches developed to this purpose. 
Some of them rely on the effects of entrainment and mixing with non-cloudy air occurring in the 
regions of the cloud near the boundaries (see, e.g., iBlythI (|l993|)). However, broad spectra ar e 
also observed inside the inner adiabatic cores, as reported by iBrenguier and ChaumatI (|200l[) . 
where no boundary effects can possibly explain their presence. Droplets themselves have been 
proposed to be sources of local variability in the environmental conditions. Indeed, droplet 
evaporation is an internal source of turbulent kineti c energy due to cooling a s sociated to the ab- 



sorption of latent heat, coupled with buoyancy (see^ Andrejczuk et aL (2004); Malinowski et al, 



(ll998h:lKorczyk et al.l (j2006)). Mor eover, the presence of droplets locally changes the water 
vapor content through phase change, IShaw et al.l (Il998h . The general conclusion is that these 
microscopic fluctuations influence the process of mixing occurring at the interface between 
cloud and clear-air. Also, the presence of ultra-giant con densation nuclei has been proposed 
to explain the large raindrops production (|Johnsonl (119821) ). Other properties such as salinity 
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and s urface curvature may produce absolute and relative broadening, as proposed in Korolevl 
(fl995 ). 

Stochastic fluctuations of the environmental conditions, induced by turbulence have been 
suggested as a broadening mechanism since the '60s when the theo ry of stochastic condensa- 
tion w as first proposed, see e.g. iMazinI (Il968|) : Levin and Sedunov ([l966t); iBartlett and JonasI 
(|1972|) . This approach explains droplet spectra broadening by observing that fluctuations of 
the ambient conditions make droplets grow at different rates. This simple idea is very pow- 
erful in that it concerns the whole cloud, regardless of the presence of additional microscopic 
mechanisms, the influence of ultra-giant nuclei and/or boundary effects such as mixing with 
dry air and entrainment. Although all these ingredients play a role in the huge natural labo- 
ratories represented by clouds, their contribution may vary according to different conditions. 
On the contrary, turbulence is a very good candidate as general broadening mechanism within 
convective warm clouds. At a Reynolds number approaching Re ~ 10^, turbulence is known to 
be highly intermittent, with statistics strongly deviating from a Gaussian one and a sub stantial 
probability of fluctuations far exceeding the standard deviation (see, e.g., iFrischI h995\) ). This 
means that droplets coming close one to the other might have previously experienced disparate 
conditions, thus invalidating any expectation based on uniform condensation. Here we show 
that some qualitative features emerging from a simple idealized setting may be well present 
in more complex models, just because they are consequences of basic properties of turbulence 
itself. In view of this purpose, under consistent hypothesis, we consider a simple but rather 
general model for condensation to be detailed below. 



Droplet growth in a fully developed turbulent flow can not be treated analytically and 
the numerical analysis becomes a fundamental tool of investigation. Numerical simulations 
of cloud turbulence are very demanding because of the huge number of active degrees of 
freedom (the microscale Reynolds number varies as a function of the spatial resolution ap- 
proximately as Rex ~ A^^^'^). As a consequence of the spatial structure of the turbulent 
fie lds, direct numerical simu lations (DNS) can either focus on large -scale structures as done 
bv ICelani et al.l dlOOSL l2007h. or resolve the small scale features as in lVaillancourt et al.l (1200 IL 
20021) : iKorolev and MazinI (|2003b : lAndrejczuk et al.l (|2004l) . The two choices are actually mu- 
tually exclusive due to finite computing resources and both have strong and weak points (see 
Section [21 for further details). Here we will try to take advantage of both these approaches, to 
achieve a deeper understanding on the role of turbulent fluctuations for the problem at hand. 
In a nutshell, we first evaluate the spreading of droplet- size distribution through a series of di- 
rect numerical simulations, at increasing resolution, matching the small scale features. Not sur- 
prisingly, each single DNS gives a small degree of s preading, as already poin t ed out in previous 
works f ocusing on small ri sing fluid parcels (see e.g. IVaillancourt et al.l(l2002|) : lKorolev and Mazin 
(l2003h : lAndrejczuk elaP (12004)). However this only tells us that turbulence at the smallest 
scales alone fails to explain the fast broadening of droplet-size distributions observed in clouds. 
Then, we evaluate the dependence of the size spectra broadening on the turbulent Reynolds 
number, i.e. on the range of spatial scales resolved. The broadening is found to increase with 
the Reynolds number of turbulence. 

Since the Reynolds numbers of real cloud turbulence are several orders of magnitude larger 
than those described by DNS, this increase must be accounted for when assessing the role of 
turbulence for condensation in clouds. In order to get more insight on this point, we analyze the 
trend of the size-spectra broadening as a function of the Reynolds number. It is quite natural 
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to derive, on the basis of consistent assumptions, expectations for the observed broadening at 
large Reynolds number. Within the model adopted, we derive upper and lower bounds on the 
trend through a dimensional analysis of the system equations. The upper bound is obtained 
by neglecting the vapor depletion due to condensation onto cloud droplets. We argue that this 
should be significant at small time lags and small Reynolds numbers, and we show that this is 
consistent with two series of DNS performed. We then reason that the lower bound, where the 
vapor fluctuations are considered immediately depleted for condensation/evaporation, should 
be more significant at high Reynolds numbers and large time lags. A further DNS matching the 
large scales of the problem suggests consistency with the argument. 

The paper is organized as follows. In Section[2lwe introduce the model for the time evolution 
of the vapor field and the droplets advected by the turbulent air flow. Section [3]is devoted to the 
numerical approach and DNS details. Results concerning the spreading of droplet-size spectrum 
are discussed in SectionHl while in Section [5] conclusions and perspectives will follow. 



2. Model equations and numerical procedures 

We focus on a turbulent velocity field advecting vapor and droplets. The latter undergo 
size changes for evaporation or condensation of the surrounding vapor. The three-dimensional 
velocity field v evolves according to the Navier-Stokes equations for an incompressible flow, 

9jv + V • Vv = -— + uAy + f , V • V = 0, (1) 

Pa 

where p is the pressure, pa is the air density and f is an external statistically homogeneous and 
isotropic forcing, providing a turbulent stationary flow. In eqlH z/ ^ 0. 15 cm^/ s is the air kine- 
matic viscosity. Note that turbulence in clouds is generated by large-scale turbulent fluctuations 
which can be strongly anisotropic. Indeed the presence of gravity - here neglected - introduces a 
preferential direction through large-scale thermal gradients and buoyancy, the engine of convec- 
tive motions. However, we can assume that for sufficiently small internal cloud cores, vertical 
stratification of the environment can be neglected and the small-scale flow is essentially forced 
by nonlinear transfer from larger scales, rather than by buoyancy. Isotropy can thus be assumed 
for these scales of motion. In Vaillancourt and Yau (2000), the authors argue that this should be 
valid in warm-cloud cores for spatial scales up to L ~ 100 m. In addition to large-scale thermal 
gradients, anisotropy is also produced through buoyancy by microscopic temperature fluctua- 
tions. Sedimenting droplets can evaporate, absorbing latent heat and thus locally cooling the 
environment. There are conditions for which t his effect can b e impo r tant for the cloud-clear ai r 



mixing at s mall turbulen t kinetic energy rates (IKorczyk et al.l (120061) : iMalinowski et al.l (|1998l) : 



Andrejczuk et al.l (|2004|) '). However, provided that we focus on the inner adiabatic core, away 
from the cloud boundaries, we can neglect this microscopic source of anisotropic fluctuations 
as a first approximation. With equation ([T]), we focus on the turbulent motion of in- cloud air, ne 



electing the role of convective motions. Note that previous two-dimensional DNS (iCelani et al 
(12005, 2007)), suggest that the qualitative effects of turbulence on condensation do not rely 
specifically on the statistical details of the turbulent regime analyzed. 

Water vapor molecules carried by the turbulent velocity field are the source for droplet 
growth by condensation. The relevant quantity for droplet condensation/evaporation is the su- 
persaturation, which quantifies the presence of vapor available for cloud particles growth. Su- 
persaturation is defined as s : = e/e^ — 1 , where e and are the vapor pressure and the saturation 
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vapor pressure respectively. Droplets are able to grow when the surrounding vapor content ex- 
ceeds the saturation point (s positive in equation moist air). On the contrary droplets tend 
to evaporate when s is negative (dry regions). Exhaustive in situ measures of the small-scale 
statistics of the vapor field are not available so far. Therefore, different models proposed cannot 
yet be validated by direct comparison with real data. 

For the sake of generality, we adopt here the simplest generalization of the well-known model 
proposed in Twomey (1959). While Twomey considered a one-dimensional equation for the 
time-dependent supersaturation function, here we consider the turbulent vapor to fluctuate both 
in space and time. For this reason we introduce the supersaturation as a field s(x, t) : = 
^(x, t) — 1, that quantifies the amount of vapor which is present at point x at time t. Since 
s(x, t) is allowed to fluctuate from positive to negative values, dry and moist regions can coex- 
ist at the same time. 

The generalization of Two mey 's equation (see alsolPrup pacher and Klett ( 1997 ): Shaw ( 2003 ): 



Korolev and MazinI (|2003[) ). for the supersaturation field is an advection-diffusion equation: 



dts + \ -Vs = kAs + Aiw - — . (2) 

Here we assume that the scalar field s(x, t) is passively advected by the turbulent flow v(x, t) 
which is not affected by its presence. In this way, we are neglecting the compositional effects 
of vapor on th e buoyancy forces acting on the flow, that are generally thought to be small (see 
StevensI ( 2005 )'). In equation (O, k ^ 10~^m^s~^ is the molecular diffusivity of water vapor in 



air and w(x, t) is the vertical component of the turbulent velocity field v. 
The term Aiw acts as a source/sink term of supersaturation resulting from the variation in 
temperature and pressure with height. It can be computed by assuming that (i) air pressure 
is quasi-hydrostatic, (ii) the pressure of dry air is similar to that of moist air, (Hi) the motion 
inside the cloud core is adiabatic (temperature changes with vertical motion according to a dry 
adiabatic lapse rate and with latent heat effects as described be low) and (iv) the atrnqspher e 



is not far from saturation . These standard arguments (see e.g. iPruppacher and Kletti (|1997|) 



Korolev and MazinI (120031) ) yield 

Lg g 



R^CpT Rq^T 



where L is the latent heat of evaporation, and Ra are the gas constants for vapor and dry 
air, Cp is the specific heat at constant pressure, and g is the gravitational acceleration. Ai can be 
interpreted as a global supersaturation gradient. 

The term —s/ts accounts for the double effect of condensation/evaporation on supersaturation: 
on one side the phase change directly modifies the water vapor content, on the other side it 
locally modifies temperature due to absorption or release of latent heat. In dry regions, s(x, t) < 
0, droplets tend to evaporate releasing vapor and cooling the environment; here, the term —s/ts 
is a source of supersaturation. Conversely in moist regions droplets tend to absorb vapor for 
condensation and to release latent heat, and — s/r^ acts as a sink term. The parameter is the 
relaxation timescale of the supersaturation and depends locally on the concentration and size of 
droplets. In each elementary volume V, Ts is defined as (see appendix A): 



n 



1=1 
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where Ri are the radii of the n droplets inside the considered volume; A2 is a function of 
thermodynamic parameters (see appendix A); is the water density; A-^ is the rate of droplet 
radius growth by condensation (see equation dH)). In the numerical procedure, we consider each 
droplet to affect the value of in the eight nodes of the grid cell surrounding its position. The 
weight of the contribution to each node is calculated via a three-linear extrapolation. Table [T] 
shows the reference values of these parameters used in the numerical experiments. 

Note that Twomey's one-dimensional equation for s is derived under a given temperature 
profile, neglecting fluctuations. It is difficult to quantify temperature fluctuations in clouds, 
since in situ microscale cloud data are unavailable. For what concerns the vapor field, this 
amounts to neglect an additional source of fluctuations induced by temperature advection and 
diffusion. The remarkable feature of the simple model identified by equations ([B and Q is 
that, despite its simplicity, it allows to identify nontrivial mechanisms leading to the spreading 
of the size spectra. 

Given the evolution equation for the Eulerian turbulent fields, we can now introduce the 
Lagrangian dynamics of cloud particles and the time evolution of their radii. A complete de- 
scription of the relation between the water vapor and the size of a droplet would imply an 
integral equation for the local dynamics occurring at the droplet surface. Note that, however, 
the typical timescales associated with the diffusional growth of an isolated droplet are much 
smaller than the faste s t time scale associated with turbulent changes in the ambient conditions 
Pruppacher and Klett (|l997|) . Therefore we can consider the droplet to be instantaneously in 



equilib rium with the surrounding vapor (a detailed quantification of this assumption is given. 



e.g., in IVaillancourt et al.l (120011) . where the authors argue that it should be valid in the con- 
densation stage, for kinetic energy dissipation rates of about e ~ 10^^ m?s^^ that we consider 
here). Under basic assumption^] we end up with the following equation for the i-th droplet 
growth rate (see also Pruppacher and Klett (1997) for further details). 

Here A3 is a function of the local conditions, air temperature and pressure, and assumed to be 
constant throughout the entire volume (variations of this parameter with temperature in typi- 
cal warm cloud conditions are smaller than 3%). According to equation dH), the growth rate 
varies from a droplet to another, since it depends on the supersaturation fluctuation s(Xj(t), t) 
measured along the trajectory Xj(t) of the single droplet. Due to turbulent transport, initially 
close droplets separate very rapidly and eventually experience disparate values of supersatura- 
tion throughout the entire cloud volume. This is the most important difference with respect to 
Twomey's model where all the droplets are exposed to the same supersaturation value. 

Cloud droplets can be described as independent. Stokes particles, whose trajectories Xj(t) 
and velocities Vj(t) evolve according to: 



'(i) Droplets are considered spherical (significant deformations from the spherical shape are typical of much 
larger drops, from sizes of hundreds of ^m); (ii) the coefficient A3 is not significantly altered by either the chemical 
composition of droplets or the size of the droplet itself (this holds for droplets already activated onto condensation 
nuclei); (iii) curvature and salinity corrections to the supersaturation at droplet surface (described by classical 
Kohler theory) a re neglected for simplicity - these are generally though t to be small for radii larger than few 



microns (see e.g. lKorolev and MazinI ( l2003h : lPruppacher and Kletti (Il997h ). 
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dt 



V,(t)-v(X,(t),t) 



(6) 



is the particle re- 



Here v(Xj(t),t) is the fluid velocity at the particle position; r^(t) — 
sponse time (or Stokes time); (3 = 3pa/ {pa + 2p^) ^ 3pa/ {'^Pw) is the air/water density ratio; 
and g is the gravity acceleration. Equations ([5]) and derived from the more general treatment 
of Maxey and Riley (1983), are valid for dilute suspensions of small spherical heavy particles. 
These hypothesis are well verified during the condensation stage, as discussed e.g. in Vaillan- 
court and Yau (2000). 

Droplet statistics and dynamics in a turbulent medium can be studied by focusing on possibly 
many different features, see e.g. Pinsky and Khain (1997) for a r eview. Recently, also thanks to 



fundamental improvements in the labora tory tracking techniques ( LaPorta et al.l(|200l|) ; lMordant et al 



(|200l|) ; lAyyalasomayajula et al.l (|2007[) ). progress in the understanding of particle motion in 
turbulence has undergone an explosive development. Together with the experimental efforts, 
theoretical analysis and numerical simulations have considered su c h different issues as par- 
ticle acceleration statistics (see e.g. lAyyalasomavaiula et al. 1/2007); 



cle spatial inhomogeneous d istribution (see e .g. lAliseda et al.i (i2002) ; Balkovsk v et al 



Be e et all (120061) ) parti- 
(p00i|); 

Collins and Keswanil i2004f) : lBec et all (120071) ). coUisional effects (see e.g. W ang et all (l2005b : 
Wilkinson et al. (200^ Granted that in the core of a warm cloud all these features can be rel- 
evant, here we focus on the particle size statistics during the early stages of the condensational 
growth. 



3. Numerical simulations and range of parameters 

As discussed earlier, DNS of cloud physics present a major problem: there is a huge number 
of degrees of freedom that cannot be described simultaneously. Turbulence is organized in 
spatial structures of typical scales ranging from the large scale L of hundreds of meters, down 
to the Kolmogorov scale r] (typically 1 mm). Similarly, the timescales range from thousands to 
fractions of a second. Within this highly turbulent medium, a population of 10^^ 10^^ droplets 
evolve. Moreover, even if droplets are much smaller than any turbulent eddy, their trajectory 
spans the whole ran ge of turbulent scales. T his yields correlations with the fluctuations of 
the vapor field - see 
as shown, e.g., in 



Celani et al] (120 05'. hooi) - and with the structures of the velocity field 
Eaton and Fessleii (il994). Therefore, turbulent motion at any scale plays a 



significant role in droplet dynamics. However, when dealing with experiments in silico, because 
of computational limitations, it is compulsory to choose a setting which describes only a limited 
range of fluctuations in the sy stem. 



Recent results, reported in lCelani et al.l (|2005Ll2007h . of direct numerical simulations in two 
dimensions pointed out the importance of the large-scale fluctuations of the vapor field. These 
provide a strongly variable environment for droplet evolution resulting in a spreading of the 
droplet-size spectrum. In such context, the small scales of turbulence cannot be resolved and 
th e analysis is limited to a s tatistically representative subset of the whole population of droplets. 
In lVaillancourt et al.l (|2002|) . the complementary setting is adopted: by concentrating on a small 
rising parcel, the authors can consistently describe the droplet evolution in full detail. This ap- 
proach provides small fluctuations that eventually produce a limited degree of spreading. Note 
that the small scales inside a cloud are the endpoint of a wide turbulent cascade involving a 
huge range of interacting spatial and temporal scales. In particular, they are not independent 
from large scales, suggesting that approaches which separate small from large scales might fail 



7 



in reproducing turbulent effects. In this respect, since DNS do not describe the whole inertial 
range of turbulence, the fluid-parcel approach does not represent a small volume inside a huge 
cloud, but rather a very small cloud. 



Here we wish to investigate how turbulence effects change with the Reynolds number by 
merging the two complementary approaches above mentioned. This is a crucial step in as- 
sessing the role of turbulence for cloud-droplet condensation, since real cloud turbulence has 
necessarily a dramatically higher Reynolds number than the simulated flows. To this purpose, 
we perform two series of direct numerical simulations at increasing resolution. The grid spacing 
Ax of each simulation corresponds to about 1 mm. By progressively increasing the number of 
grid points for each spatial direction, we can resolve larger integral scales L, defining the size 
of the cloud. More precisely, we consider two series of four numerical experiments, labeled as 
run (a), (b), (c) and (d), with 64^, 128'^, 256'^ and 512'^ grid points respectively. The two se- 
ries of simulations have different initial liquid water content (LWC): ~ 1.2 g/cm^ for series 1 
and ^ 0.07 g/cm^ for series 2. The integral scale of the system va ries from L ^ 9 cm up to 



L ~ 70 cm. The microscale Reynolds numbers Rcx ^ Vl5 Re (see FrischI (|l995|) ) range from 



Re\ ~ 40 to Rcx ~ 185, which represents the state-of-the-art for DNS in cloud physics. 
The ratio between the air kinematic viscosity and the vapor molecular diffusivity, also called the 
Schmidt number, is Sc = u/k = 1, so that the flow and the scalar dissipative scales are of the 
same order. Table [2] lists all the relevant DNS parameters. Clearly the process of doubling the 
resolution, if iterated, would ideally lead to the description of the whole range of scales from 
1] ^ 1 mm to L ^ 100 m. As we can only perform the first few iterations of this process, the 
next step will be to discuss if this can be extrapolated to give information on Reynolds numbers 
which are not accessible by DNS. An attempt in this direction will be discussed at the end of 
next section with some detail. 

For each DN S, equat i ons ([B and (O are integrated using pseudospectral methods with 2/3- 



rule de-aliasing (lOrszagI (|197 1)). in a three-periodic box. Molecular viscosity (and diffusivity) 



is chosen so as to match the Kolmogorov lengthscale with the grid spacing t] ~ 6x: this choice 
ensures a good resolution of the small-scale dynamics. Kinetic energy is injected at an av- 
erage rate e , by k eeping constant the total energy in each of the first two wavenumber shells 



(|Chen et al.l (119931) ). The scalar field, also integrated on a triply periodic box, is forced by the 



assigned gradient Ai; the term — s/r^ does not contribute as long as droplets are not injected 
into the flow. Time stepping is done using a 2"*^ order Adam-Bashfort scheme and the time step 
is chosen to accurately resolve the smallest turbulent fluctuations and the particle acceleration. 
We obtain a statistically stationary state for the velocity and supers aturation fields with no 
droplets by integrating equations ([T]) and Q for few large-scale eddy turnover times = 
L/vrrns- Figurc [T] shows the supersaturation spectrum at the stationary state for run (d), Rcx ~ 
200, before p articles injection. In agreemen t with classical Kolmogorov-Obukhov-Corrsin the- 



ory (see e.g. iTennekes and Lumleyl (Il972|) ). this exhibits a k^^^^ power law behavior in the 
Fourier space. Since the scalar spectrum is peaked on the large scales, as the integral scale 
increases, we approach larger and larger fluctuations. We indicate with the supersaturation 
standard deviation in the stationary state before droplet injection. In the inset of Figure [B we 
show that (T° increases linearly with the size of the system as expected from a dimensional bal- 
ance of terms in equation yielding cr° ~ AiL. Such increase can be understood from the 
physical viewpoint: larger cloud sizes correspond to larger displacements and stronger adiabatic 
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cooling. This directly provides larger fluctuations in the vapor field through the term Aiw in 
equation (|2l). Without the feedback coming from droplets' absorption of vapor, the value cr° 
would remain as in the stationary state, i.e. ~ AiL ~ AiWrmsTh- 

Once the steady state has been attained, a monodispersed population of droplets (with initial 
radius Ri = 13 fxin for series 1 and Ri = 5 fim for series 2) is injected into the flow. Droplet 
concentration is for all experiments ^ 130 cm^^, which means that, on average, about one out 
of 2 ^ 4 cells contains a droplet. In the largest simulation we followed the time evolution of 32 
millions droplets. Initially, these are distributed randomly in space, according to a statistically 
homogeneous Poisson distribution. Lagrangian equations ([5]) and © for particle motion are 
integrated simultaneously with those for Eulerian fields © and ©. Particle initial velocities 
are set equal to the local fluid velocity H. To obtain droplet velocity from equation Q, the 
underlying flow velocity at the particle pos ition has to be compute d. This is done via a linear 
interpolation in the three spatial directions (|Yeung and Popd (|i988|)), which was demonstrated 
to be adequate to obtain well-resolved particle acceleration. Similarly, we compute the vapor 
field s(X(t)) at the particle position. 

Coupling between droplets and the vapor field takes place via the term — s/r^ of eq. ©, so 
that the initial vapor available is consumed for condensation onto droplet surface or released 
for droplet evaporation. In typical cloud conditions, is of order 1^10 s, much smaller than 
typical large-scale eddy turnover times. As a result, the initial supersaturation standard deviation 
is depleted from cr° ~ AiWrmsTi to a lower value, that can be estimated by a dimensional 
analysis of equation ~ AiWrms'^s- We pointed out above that grows with the 

size of the cloud and the Reynolds number: note that a"^ grows with Reynolds as well, since 
we can estimate Wrms ~ Vr^Re^"^ , where is the velocity fluctuation at the viscous scale. 
Considering that ~ alTg/Ti, the ratio between the large-scale eddy turnover time and 
the supersaturation absorption time is important to evaluate the supersaturation fluctuations 
available for condensation of cloud droplets. The complete system -flow, scalar and droplets- 
has been studied, at the largest resolution, for about two large-scale eddy turnover times. Longer 
time integrations were performed at lower resolutions. 



4. Results and Discussion 

Our simulations start with a spatially uniform Poissonian distribution of monodispersed 
droplets of radius Rq = 13 /im for series 1 and Rq = 5 fim for series 2, and with vanishing 
acceleration. As the cloud particles are released, they explore the entire volume and experience 
the range of vapor fluctuations available in the system, of initial standard deviation cr°. During 
an initial transient the fluctuations of the vapor field decrease due to the feedback of droplets. 
In Figure [2l the droplet square size distributions P{R^) are shown for runs (a)-(d) after one 
large-scale eddy turnover time for series 1 and series 2. A small degree of spreading is present 
for each simulation, increasing with the size of the cloud % This is due to the fact that - as 
already discussed above - larger domain sizes correspond to larger supersaturation fluctuations 
and longer large-scale eddy turnover times. In other words, when evolving in a larger cloud 
volume, droplets are exposed to more and more intense fluctuations of vapor, for longer and 

^Tests have been performed by setting the droplet initial velocity equal to the terminal velocity vt = gra- We 
did not observe any significant deviations in the results. Indeed particles rapidly equilibrate to the flow on a time 
scale of the order of their response time. 

^Note that we impose a vanishing mean value for the initial supersaturation field, thus a constant mean droplet 
radius. 
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longer times. The knowledge of these characteristic times turns out to be crucial, as discussed 
in the sequel. 

For t < Tl, the standard deviation a/ja of the square radii -R^(t) increases linearly in time as 
shown in the inset of Figure [21 This means that at short time lags droplet surface grows with the 
vapor fluctuation initially experienced and does not feel the underlying local variations 0. After 
few absorption times have elapsed, the vapor field standard deviation is depleted from cr° to 
(j^, thus slowing down the broadening of droplet size distribution. The mean vapor absorption 
time, estimated through equation (|3]) with average concentration ^130 drops cm^^ and mean 
radius 13 /im and 5 fxm, turns out to be ~ 2.5 s for series 1 and ^ 7 s for series 2. Since this 
time scale is comparable with the eddy turnover time for simulations (a)-(d) (see table [T] and |2l), 
the supersaturation fluctuations do not change considerably from their initial value during T^. 
For this reason, the linear growth of 0-^2 in time extends to the whole large-scale eddy turnover 
time Tl and the final spreading, measured att = Tl, is well approximated by: 

a^2(Tz.) ~2A3a°Ti.. (7) 

The square radius standard deviation aLi'2{TL) is found to be similar for the two series of sim- 
ulations, pointing once more to a limited role of vapor absorption by condensation when the 
Reynolds number is small and ~ T^. 

In Figure [3l the standard deviation o"i?2(Tj^) is shown as a function of the microscale Reynolds 
number Rex characterizing flows (a), (b), (c) and (d). We immediately observe that the final 
broadening grows with the Reynolds number for both series of simulations. Given the following 

dimensional relations Tl ~ TrjRex, Wrms ~ Vr]Re]!'^ and cr° ~ AiWrmsTL, it is straightforward 
to translate the simple expectation (|7]) in a scaling relation with the Reynolds number: 

aR2 ~ AsA^v^T^Re'/' , (8) 

based on self-similarity of the growth process during the turbulent regime. In Figure [3] we show 
for comparison the dimensional expectation ([8]). 

It is natural to wonder whether this trend can give information on the final broadening 
achieved at Reynolds numbers higher than Rex = 200. This delicate point deserves insight 
since the Reynolds number of real cloud turbulence typically achieve huge values that cannot 
be described by direct numerical simulations. 

At larger Reynolds number and large-scale eddy turnover times, droplets significantly absorb 
the surrounding vapor, so that o"° is no longer a good approximation of cr^. This means that 
the scaling ^ is an upper bound for the actual spectral broadening at larger Reynolds number. 
A lower bound can be simply obtained by replacing the supersaturation fluctuations with their 
long time value cr^ ~ a'^Tg/TL, yielding a correction to the upper bound ([8]) of a factor Ts/Tl- 
Note that the scaling relation in ([8]) would also change into a slower growth of the final broad- 
ening with the Reynolds number, proportional to Re^x^. 

The ratio Ts/Tl fluctuates both in space and time according to the local properties of both 
turbulence and droplet population. Given the Kolmogorov timescale, we can estimate the 
large-scale eddy turnover time at larger Reynolds numbers through Tl ~ r^^Rex, as men- 
tioned above. For a cloud of typical size L ~ 100 m, liquid water content LWC 1.2 g/m? with 

"^A proportional linear growth is also observed for the standard deviation of droplet radii R{t), (not shown). 
This is because, given the tiny supersaturation fluctuations, the droplet-size distribution is close to a Gaussian and 
the mean radius is much larger than the standard deviation. Hence cr^a ~ 2ct^ + 4cr|j.(i?)^ « Aa\{B)^. 



10 



130 drops cm~^ and Reynolds numbei|f| Rex ~ (4000 ^ 7000), the large-scale eddy turnover 
time is ~ 150 s while does not depend on Reynolds and its mean value is ^ 2.5 s, 
as in series 1. The extrapolation of the upper bound ([8]) together with the lower bound gives: 
(3.3 ± 1.6)i2m'^ < aR2{TL) < (200 ± im)nm^. 

Validating these predictions by a direct numerical simulation, describing the proper number 
of droplets and the whole range of space-time scales involved, is of course not possible. Instead, 
we present a direct numerical simulation, labeled run (e), designed so as to reproduce the large- 
scale parameters typical of a cloud of size L = 100 m: Tl ^ 150 s, cr° ^ 2%, Vrms ~ 0.6 ms^\ 
L fii 100 m. The goal of this further experiment is to verify that at long space and time scales 
(Tl is now much larger than Ts), the lower bounds on vapor fluctuations ag and on the square 
radius fluctuations aR2 become more relevant. Of course, the small-scale parameters do not 
match the realistic ones: in this run the smallest resolved scales are i] ^ 25 cm and ~ 4 s, 
and the turbulence Reynolds number is Rex = 185. The initial radius of droplets is 13 fim and 
the numerical resolution is N'^ = 256^ grid points. Space-time integration of the system has the 
same features described in Section [3l and reference values for the physical parameters entering 
the model equations are those listed in Table 1 . 

Clearly, we can not follow the evolution of 130 cm^^ droplets, since in a volume of (100)^ 
they would sum up to A^* = 1 .3 x 10^"*. The traditional cloud physics approach avoids this prob- 
lem by focusing on a nun iber density function rep r esenting the local co ncentration of droplets 



with a given size (see e.g. lAndrejczuk et al.l (|2004|) : lJeffery et al.l (120071) ). However, a complete 



continuous descripti on for the turbulent transport of inertial particles is still an open issue (see 



Boffetta et al.l (|2007[) and references therein). Though more computationally demanding, a La- 
grangian approach turns out to be more appropriate. We choose then to consider the complete 
evolution of a subset N^rops of droplets, representative of the whole population. Of course, 
the whole population would absorb for condensation more vapor than the representative subset 
does. Therefore an algorithm is needed to account for the correct feedback of droplets on the 
supersaturation field. In order to accomplish this task, we simply normalize the field with 
a factor N* /Ndrops- Details on the meaning and the convergence of this simple algorithm are 
discussed in appendix B. 

The supersaturation fluctuations for simulation (e) start with ~ AiWrmsTi ~ 2%; after 
few vapor absorption times they are depleted to af" ~ 0.04% and oscillate around this value 
for all later times. The relevant spreading of the square size distribution shown in Figure |4] can 
be quantified in terms of the standard deviation of the radius and of the square radius. After 
one large-scale eddy turnover time, we obtain aRiT^) — (0.30 ± 0.04) /xm and or2{Tl) ^ 
(7 ± 1) yum^, respectively. Although simulation (e) cannot resolve the whole range of spatial 
scales of turbulence, by matching the parameters on the large scales of the problem it is able to 
reproduce the intensity of the large-scale fluctuations. 

The final spreading achieved is very close to the lower bound, obtained above by assuming 
that the supersaturation fluctuations stabilize to af and rapidly forget their initial condition. 
This is consistent with the picture drawn and supports the observation that the lower bound is 
more relevant when dealing with large space-time scales. The result suggests that this could 
be the case for large Reynolds number turbulence. Let us remark however that care must be 
taken in this respect, since the spatio-temporal complexity of turbulence at Rex ~ 5000 is not 
described by the direct numerical simulation (e) presented. 



^To put errorbars on the Reynolds number we consider fluctuations of 25% on the average value, the same level 
of fluctuation observed in simulation runs (a)-(d). 
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5. Conclusion and Perspectives 

Turbulent fluctuations have been shown to play a role for the broadening of the droplet- size 
distribution in an idealized setting for condensation in warm adiabatic cloud cores. A qual- 
itative explanation for this observation relies on the non-trivial spatial structure of turbulent 
fields where droplets collect strongly different histories. The droplet size spectra broadening is 
shown to depend on the actual range of spatial scales characterizing the turbulent structure of 
the fields. In particular, we show an increase of broadening with the microscale Reynolds num- 
ber by means of two series of direct numerical simulations up to Rex ~ 200 and with millions 
of droplets. 

The general outcome of this result is that the role of turbulence for condensation in clouds may 
not be assessed just by focusing on small Reynolds numbers, where a limited range of spatial 
scales is considered. Therefore, a strategy must be conceived to extrapolate the results of DNS 
or laboratory experiments, that cannot achieve the huge Reynolds numbers of real cloud turbu- 
lence. This is also the case for fluid parcel models which in fact do not represent a portion of a 
large cloud, but rather a very small cloud. Indeed, it is well-known that initially close droplets 
inside a parcel would not remain close for a long time, but separate explosively due to turbulent 
transport. 

In particular, though at the Reynolds numbers considered the broadening is relatively small, it 
may potentially achieve relevant values at the huge Reynolds numbers of real cloud turbulence. 
In order to get a better intuition on this potential, we perform a dimensional analysis of the final 
broadening viewed as a function of the Reynolds number. At moderate Reynolds numbers, it 
is possible to neglect the feedback of droplets on the vapor field. This yields a dimensional 
expectation consistent with the DNS data. 

At larger Reynolds numbers, however, droplets are expected to consistently modify the sur- 
rounding vapor field. Therefore, the extrapolation of the dimensional behavior provides an 
upper bound to the actual spectral broadening at large Reynolds numbers. A correction is pro- 
posed to estimate the effective absorption of vapor for condensation at large Reynolds numbers, 
yielding a lower bound on the spectral broadening. A further simulation designed to discuss 
the validity of this argument shows that at large space-time scales the behavior is well approx- 
imated by the lower bound. This suggests, within the limits of simple dimensional arguments, 
that the lower bound could be a good candidate to estimate the size-spectra broadening at large 
Reynolds numbers. 

In this work, we focused on the role of turbulence for droplet condensation at moderate- 
high Reynolds numbers. In order to enlighten the basic mechanisms yielded by turbulence, 
we considered a simple model describing the essential physical processes. More sophisticated 
models have been proposed in the literature accounting for additional microphysical and ther- 
modynamic couplings. Ingredients such as the explicit dependence of the vapor field on the 
temperature fluctuations, the microscopic interactions of droplets with the turbulent fields and 
the role of buoyancy both at small and large scales may yield corrections to the presented re- 
sults. Other aspects, such as droplet preferential concentration and modified relative velocity 
effects, were not specifically addressed in this work, though they are in fact included in the 
model. These have been associated to th e efficiency of collision/coalescence processes (see 



Wang et al.l (|2005l) : [Wilkinson et al.1 (120061) ). before they become dominated by gravity. 
The problem here considered clearly deserves further theoretical, experimental and numerical 
insight and represents a promising challenge for future research. In particular, the ideas at the 
basis of this work could be well explored through laboratory experiments. The higher Reynolds 
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numbers covered would allow to verify and complement the picture drawn through the numeri- 
cal analysis presented. 
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APPENDIX A 

Water- vapor interaction 

The expression of the absorption time can be computed directly from the classical form 
of Twomey's model by identifying the terms: 

S . dpL . RaT SwaC? 

— = ^2-7-, where A2 = h — — (Al) 

where Pl = m^/V is the density of liquid water of mass inside volume V; Syja is the ratio 
between the molecular weight of water and dry air; C is the latent heat of water evaporation; 
Ra is the specific gas constant for dry air; T is the absolute temperature and Cpa is the specific 
heat of dry air at constant pressure p (see Pruppacher and Klett (1997)). The evaluation of the 
constant A2 in typical warm cloud conditions gives the reference value shown in Table [Tl Its 
fluctuations with temperature are less then 1%, so that we assume A2 to be constant. 
Since we focus on inner cloud cores, we neglect entrainment; then droplets inside the volume 
V are the only responsible for the local change in vapor content: 

dpL 1 dmyj 1 " ^dRi Anp^As ^ 

- T7Z^47rp^i?i— = — - — }_^RiS, (A2) 



dt V dt V fr[ " dt V 



i=l 



where (see eq. (13-28) of Pruppacher and Klett (1997), ignoring small curvature and solute 
terms for activated droplets), 

^3=f^^+ 2 ) (A3) 

where R^ is the specific gas constant for moist air, is the density of water, Dy is the thermal 
diffusivity of water vapor in air, and ka is the air thermal conductivity. 

In eq. (IA2I) . i?j are the radii of the n droplets inside the volume V. The rate of variation of the 
radius is given by equation (U), where s is considered equal for each droplet inside the small 
volume V. In the numerical analysis, is a cube of edge r], similar to the vapor diffusive scale. 
The fluctuations of the scalar field are tiny under this scale and all the droplets inside the volume 
V experience approximately the same value of s. In this sense, the Twomey's parcel, where no 
spatial fluctuations of the supersaturation field are accounted for, corresponds to our grid cell. 
From equations (lAll) and (|A2I) we end up with the expression (|3]) for the absorption time r^. 
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APPENDIX B 



Renormalization of droplets population 

Simulation (e) is designed to describe a volume of (100 m)'^ with 130 cm~^ droplets. A 
major drawback of this simulation is that it can not describe turbulence structure at the smallest 
scales, below i] = 25 cm and = 4 s. Also it can not follow the individual history of the 
correct number of droplets. This might influence the statistics of vapor fluctuations since the 
correct feedback of droplets on the vapor field is not accounted for. In particular, a concen- 
tration of 130 cm~^ droplets yields A^* ~ lO^'' total number of particles, whereas we consider 
the evolution of several millions droplets (Ndrops) due to computational limits. Clearly, 10^^ 
droplets would absorb much more vapor for condensation than a few millions do, thus an algo- 
rithm is needed in order to account for this. The simplest way to estimate the feedback of the 
whole population on the supersaturation field is to normalize with a factor Ndrops/N*. This 
amounts to consider each droplet as representative of N* /N drops equal particles in the same 
volume {5xY, where 5x ^ 25 cm. 

To evaluate the reliability of this approach, we must consider the single grid cell. If there is 1 
representative drop per cell, the algorithm computes the feedback on vapor as if the average size 
of droplets in the cell was exactly the radius of that unique representative. On the contrary, a 
volume of (25 cm)^ - the grid cell size in simulation (e)-, should contain several droplets span- 
ning a whole size spectrum whose mean value is not well represented by one single droplet. Of 
course, a higher number of particles would better represent the local mean radius and we expect 
the algorithm to converge for density values larger than 1 droplet per cell. We test this expec- 
tation, by repeating simulation (e) for a time lag t = T^/A, with 0.4, 1, 3 and 6 droplets per 
grid cell, which correspond to a total number of Ndrops ~ 7, 17, 50 and 100 millions droplets, 
respectively. Consistently, we use different renormalization factors N* /N drops- On average, the 
supersaturation absorption timescale is the same in the four simulations (not shown), pointing 
to a correct renormalization. 

In figure [51 left panel, the four standard deviations of the supersaturation field are shown 
in time. They all start with the same initial value (about AiWrmsTi ~ 2% - not shown), and 
relax to a value of the order of AiWrms'Ts, with some differences for the different concentra- 
tions of representative droplets. These deviations decrease when the number of representative 
droplets increases, showing that the algorithm renormalizing the feedback of droplets on va- 
por converges. Consistently, the evolution of the size spectra broadening (shown in figure [5J 
right panel) tends to collapse to a unique curve when the number of droplets per cell is larger 
than 1. According to this analysis, we choose to perform the complete simulation lasting one 
large-scale eddy turnover time with 1 and 3 droplets per cell, which correspond to a total num- 
ber of about 17 and 50 millions droplets, respectively. The error bar on the final broadening is 
estimated from these two runs. 
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List of Figures 



Log-log plot of the stationary supers aturation spectrum for the run (d) at Re\ ~ 
200. It shows a k"^/^ power law behavior, as expected from Kolmogorov- 
Obukhov-Corssin theory. The turbulent velocity field also displays a Kolmogorov 
spectrum (not shown). Inset: log-log plot of the standard deviation of the super- 
saturation field (7°, measured in the stationary state, versus the size of the system 
L. The behavior is in agreement with the dimensional prediction cr° ~ AiL. . . |20 
Droplet square size distribution P{R^) measured after one large-scale eddy 
turnover time T^, for the 4 DNS of series 1 (left panel) and series 2 (right 
panel). At increasing the turbulence Reynolds number, we move from inner to 
outer curves. Symbols: run (a) red triangles; run (b) green dots; run (c) blue 
circles and run (d) purple squares. In each run, droplets' initial size distribu- 
tion (not shown) h 5{R — Rq), with _Ro = 13 iim for series 1 and Rq = 5 fxm 
for series 2. Each simulation presents a small degree of spreading, which in- 
creases with the Reynolds number. Inset: log-log plot of the time evolution of 
the standard deviation of the square size distribution aji2 for the runs (a)-(d) 



(from bottom to top). Symbols are the same of the main frame |21 

Log-log plot of the spreading of droplet size distribution aji2 [Ti] for the square 
radius R^, measured after one large-scale eddy turnover time T^, as a function 
of the Reynolds number Rex- Data refer to simulations (a)-(d) of series 1 (red 
squares) and series 2 (green circles). For both series, the spreading is larger 
as the Reynolds number increases, since droplets evolve in conditions which 
are more and more differentiated. The dimensional prediction app. ~ ciRe\ 
with ^ = 5/2, ci oc AiA^VriT^^ is shown for comparison (see the text). The 
extrapolation of this law gives an upper bound to the final spreading for the 
target cloud of parameters L = 100 m, a" ~ 2 % and Rex ~ 4000 ^ 7000: 
(7^2* ~ (200 ± 100)/im^. A correction of the expectation accounting for vapor 
depletion due to droplet feedback gives a^l* ~ (3.3 ± 1.6) fim'^ for series 1, and 

(7^2* ~ (9.3 ± 4.5) yum^ for series 2 (see the text) |22 

Droplet square size distribution P(i?^) measured after one large-scale eddy 
turnover time ~ 150 s for run (e), matching the large scale cloud parameters. 
In the inset, time evolution of the standard deviation of the radius distribution, 
aR{t), and of the square radius distribution, (Tr2 (t). At time t = Tl, we measure 



aniTL) ~ (0.30 ± 0.04) /im and (Tr2{Tl) ~ (7 ± 1) /im^ |23 

Left panel: Supersaturation standard deviation as a function of time. Right 
panel: Droplet-square-radius standard deviation crj^2 as a function of time. Re- 
sults are shown for simulation (e), with 4 different numbers of droplets N drops 
representative of a population of A^* ^ lO^'' droplets. The feedback of the 
whole population on the supersaturation field is accounted for by renormaliz- 
ing Tg with a factor N* /N drops (see text). The algorithm converges when the 



number of representative droplets is larger than 1 / cell |24 
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Fig. 1. Log-log plot of the stationary supersaturation spectrum for the run (d) at Rex ~ 200. 
It shows a k~^^^ power law behavior, as expected from Kolmogorov-Obukhov-Corssin theory. 
The turbulent velocity field also displays a Kolmogorov spectrum (not shown). Inset: log-log 
plot of the standard deviation of the supersaturation field cr°, measured in the stationary state, 
versus the size of the system L. The behavior is in agreement with the dimensional prediction 
a° ~ All. 
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Fig. 2. Droplet square size distribution P{R^) measured after one large-scale eddy turnover 
time Tl, for the 4 DNS of series 1 (left panel) and series 2 (right panel). At increasing the tur- 
bulence Reynolds number, we move from inner to outer curves. Symbols: run (a) red triangles; 
run (b) green dots; run (c) blue circles and run (d) purple squares. In each run, droplets' initial 
size distribution (not shown) is 5{R — Rq), with _Ro = 13 fim for series 1 and Rq = 5 /im for se- 
ries 2. Each simulation presents a small degree of spreading, which increases with the Reynolds 
number. Inset: log-log plot of the time evolution of the standard deviation of the square size 
distribution aji2 for the runs (a)-(d) (from bottom to top). Symbols are the same of the main 
frame. 
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Fig. 3. Log-log plot of the spreading of droplet size distribution (7^2 (T^) for the square radius 
R"^, measured after one large-scale eddy turnover time T^, as a function of the Reynolds number 
Rex. Data refer to simulations (a)-(d) of series 1 (red squares) and series 2 (green circles). For 
both series, the spreading is larger as the Reynolds number increases, since droplets evolve in 
conditions which are more and more differentiated. The dimensional prediction aR2 ~ Ci-Re^ 
with ^ = 5/2, ci oc AiA^VnT^ is shown for comparison (see the text). The extrapolation of this 
law gives an upper bound to the final spreading for the target cloud of parameters L = 100 m, 
a° ~ 2 % and Rex ~ 4000 ^ 7000: a^f ~ (200 ± 100)fim^. A correction of the expectation 
accounting for vapor depletion due to droplet feedback gives ~ (3. 3 ±1.6) fim? for series 1, 
and cr^2* ~ (9.3 ± 4.5) /im^ for series 2 (see the text). 
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Fig. 4. Droplet square size distribution P{R^) measured after one large-scale eddy turnover 
time Ti ~ 150 s for run (e), matching the large scale cloud parameters. In the inset, time 
evolution of the standard deviation of the radius distribution, crji(t), and of the square radius 
distribution, aR2{t). At time t = Ti, we measure <yR{Ti) ~ (0.30 ± 0.04) jim and aR2{TL) ~ 
(7± l)/im2. 
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Fig. 5. Left panel: Supersaturation standard deviation as a function of time. Right panel: 
Droplet-square-radius standard deviation 0^2 as a function of time. Results are shown for 
simulation (e), with 4 different numbers of droplets Ndrops representative of a population of 
N* ^ 10^^ droplets. The feedback of the whole population on the supersaturation field is ac- 
counted for by renormalizing with a factor N* /N^rops (see text). The algorithm converges 
when the number of representative droplets is larger than 1/cell. 
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List of Tables 



1 Reference values for the physical parameters used in the numerical experiments . 
Ai is the global supersaturation gradient. A2 and are functions of the am- 
bient thermodinamic parameters (see appendix A): as reference values for the 
temperature and pressure we used T = 283/^ and p ^ lOOOhPa. The values 
of the vapor relaxation time r^, the droplet radius R and the Stokes number St 
correspond to averages on the initial condition of droplet population 

2 Parameter of the DNS, series 1 and 2. From left to right: number of gridpoints 
A^^, integral scale L, large-scale eddy turnover time T^, microscale Reynolds 
number Rex, average kinetic energy dissipation rate e, Kolmogorov spatial scale 
7], Kolmogorov timescale r^, initial supersaturation standard deviation cr^, ve- 
locity standard deviation Vrms and number of droplets N drops 
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Table 1. Reference values for the physical parameters used in the numerical experiments. Ai 
is the global supersaturation gradient. A2 and A3 are functions of the ambient thermodinamic 
parameters (see appendix A): as reference values for the temperature and pressure we used 
T = 283K and p ^ lOOOhPa. The values of the vapor relaxation time r^, the droplet radius R 
and the Stokes number St correspond to averages on the initial condition of droplet population. 



label 


Ai 


A2 


A3 
{/irn^s'^) 


drops / ^ 

(cm~^) 


LWC 
(g/m^) 


is) 


R 


St 


series 1 


5 X 10-^ 


350 


50 


130 


1.2 


2.5 


13 


3.5 X 10-2 


series 2 


5 X 10"^ 


350 


50 


130 


0.07 


7 


5 


5 X 10"^ 
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Table 2. Parameter of the DNS, series 1 and 2. From left to right: number of gridpoints N^, 
integral scale L, large-scale eddy turnover time Tl, microscale Reynolds number Rex, average 
kinetic energy dissipation rate e, Kolmogorov spatial scale 77, Kolmogorov timescale r^, initial 
supersaturation standard deviation (7°, velocity standard deviation Vrms and number of droplets 

^drops • 



label 




L 


Tl 


Rex 


e 


V 












drops 






{cm) 


is) 






{cm) 


{s) 


{%) 




{ms~^ 


) 


(xlO^) 


(a) 


64^ 


9 


2.0 


40 


10^3 


0.1 


0.1 


1.5 X 10^ 


-3 


4 X 10" 


-2 


0.93 


(b) 


128^ 


18 


3.5 


65 


9. X 10"^ 


0.1 


0.1 


3.4 X 10- 


-3 


5. X 10- 


-2 


8.2 


(c) 


256^ 


38 


5.5 


105 


10-3 


0.1 


0.1 


6.1 X 10" 


-3 


7. X 10- 


-2 


71.2 


(d) 


512^ 


70 


7.6 


185 


1.1 X 10-3 


0.1 


0.1 


1.2 X 10- 


-2 


1. X 10- 


-1 


320 
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